Physics of Synchronized Neutrino Oscillations Caused by Self-Interactions 
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In the early universe or in some regions of supernovae, the neutrino refractive index is dominated by 
the neutrinos themselves. Several previous studies have found numerically that these self-interactions 
have the effect of coupling different neutrino modes in such a way as to synchronize the flavor 
oscillations which otherwise would depend on the energy of a given mode. We provide a simple 
explanation for this baffling phenomenon in analogy to a system of magnetic dipoles which are 
coupled by their self-interactions to form one large magnetic dipole which then precesses coherently 
in a weak external magnetic field. In this picture the synchronized neutrino oscillations are perfectly 
analogous to the weak-field Zeeman effect in atoms. 

PACS number(s): 14.60.Pq, 



I. INTRODUCTION 

In a two-flavor system of mixed neutrinos, the flavor 
content of a given state osciflates with the frequency 



Am /2p where Am 



ml 



is the neutrino mass- 



squared difference and p is the momentum. Therefore, if 
a neutrino ensemble encompasses many modes with many 
different momenta, these modes develop growing relative 
phases so that the overall flavor content of the ensem- 
ble quickly decoheres. This trivial effect is illustrated in 
Fig. 1^ for an ensemble of neutrinos (no anti-neutrinos) 
with a thermal momentum distribution at temperature 
T. The vacuum mixing angle was taken to be sin 29 = 0.8 
and all neutrinos were originally in a pure Ve state. In 
our example the momentum distribution is very broad so 
that the flavor decoherence takes place within about one 
oscillation period (dotted line in Fig. |^). 

This behavior changes dramatically when the neutri- 
nos feel a signiflcant weak-interaction potential caused 
by the presence of the other neutrinos. We express the 
strength of the neutrino-neutrino potential in terms of 
the parameter 




2\/2Gf?t-i/Po 



(1) 



where po = (p~^}~^- When the neutrino-neutrino poten- 
tial is comparable or much larger than a typical Am^/2p, 
corresponding to k = C(l) or larger, the modes get 
locked to each other — the entire ensemble oscillates with 
a common frequency Wsynch which corresponds to a cer- 
tain average of Am^/2p (Fig. |l|). This stunning effect 
was first discovered in numerical studies of early-universe 
neutrino oscillations ^ and then elaborated and applied 
in a large series of papers |||-^ . 

We note that the mode synchronization effect discussed 
in pfl] is unrelated to our present case. When frequent 



FIG. 1. Total Ue survival probability as a function of time, 
where r = (Am^/2po)t and po = {p~^)~^ ^ 2.2 T. The 
curves are for different values k of the neutrino self-coupling 
as indicated where k = corresponds to vacuum oscillations. 



flavor-blind collision occur on a time scale much faster 
than the oscillation period, then the energy of a given 
neutrino is averaged over an oscillation period, leading 
to a common oscillation frequency. 

Our present mode synchronization effect illustrated in 
Fig. 0is a strictly nonlinear effect caused by the neutrino- 
neutrino self-interactions and as such seems difficult to 
understand. We presently develop a very simple and 
physically transparent theory of this effect, taking full ad- 
vantage of the equivalence of our problem with the spin- 
precession of a magnetic dipole in magnetic fields. It will 
become clear that the essential effect of the neutrino self- 
coupling is to lock the individual neutrino modes to form 
one large "magnetic moment" which then spin-precesses 
in a weak external field. Therefore, the equation of mo- 
tion returns to a simple linear form. 

Our approach provides a transparent and intuitive an- 
alytic framework which nicely accommodates and illumi- 
nates the results of the previous literature which were 
largely based on numerical studies. Once the framework 
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is established, it is easy to study various generalizations 
and special cases that otherwise would be difficult to pre- 
dict or understand. 

In Sec. II we begin by setting up the equations of 
motion for neutrino oscillations with the inclusion of a 
neutrino-neutrino potential. In Sec. Ill we develop the 
picture of coupled magnetic moments to explain the syn- 
chronization effect. In Sec. IV we include anti- neutrinos, 
a situation where the system can behave qualitatively dif- 
ferent from the neutrino-only case. Finally we summarize 
our findings in Sec. V. 



II. EQUATIONS OF MOTION 

Our starting point is the well-known spin-precession 
picture for neutrino oscillations in vacuum |[Tl|-[l^ 



2p 



B X P. 



(2) 



Here, P is the polarization vector in flavor space of a 
neutrino mode with momentum p. In the usual way, the 
z-component of P gives us the probability for finding 
the neutrino, say, in the electron flavor state by virtue of 
prob(;/e) = ThevectorB = (sin 26*, 0, - cos 26*) 

with the mixing angle 9 gives us an effective "magnetic 
field" around which P precesses. Therefore, P plays 
the role of an angular momentum vector while M = 
(Ato^/2p) P plays the role of a magnetic dipole moment 
associated with P. The quantity ATn'^/2p = |M|/|P|, 
gives us the proportionality between M and P and thus 
plays the role of the "gyromagnetic ratio" for a given 
mode, determining the rate of precession. 

In much of the literature, the equation of motion is 
written in the form P = V x P without distinguish- 
ing clearly between the effective angular momentum P 
and its associated magnetic moment. For our present 
discussion this distinction is crucial. Still, we could 
split V in different ways between the gyromagnetic ra- 
tio, the unit of magnetic moment and the B-field 
which really stands for fiB. For example, we might 
have used po /p as the gyromagnetic ratio with po some 
typical or average momentum, and defined B — /iB — 
(AtoV2po) (sin 26*, 0, - cos 26*). However, we have pre- 
ferred to avoid introducing an additional quantity pq, and 
it is convenient to define B as a unit vector. 

We will frequently consider the polarization vector for 
an entire ensemble of neutrinos. 



(3) 



i.e. we consider a large volume V filled homogeneously 
with Ni, neutrinos. We also assume an isotropic distri- 
bution of the momenta so that it suffices to specify the 
modulus of the momentum of a given mode, pj = \pj\. 



There is no closed equation of motion for J because the 
individual modes oscillate with different frequencies. Ev- 
idently, however, the projection of J on B is conserved. 
On the other hand, the fast precession of the individ- 
ual mode polarizations around J average the transverse 
components of the individual modes to zero so that the 
asymptotic value Joo = (B • J)B obtains. (Recall that in 
our definition B is a unit vector.) For maximum mixing, 
and if initially all neutrinos were in a fiavor eigenstate, 
Joo — 0, corresponding to an incoherent mixture of both 
flavors. 

A background medium consisting, say, of protons, neu- 
trons and electrons modifies the "magnetic field." As- 
suming our two-fiavor system involves electron neutrinos 
the substitution is 



A?- 



2p 



B 



At 



2p 



■ B + V2G 



F rig z 



(4) 



with Tie the electron number density. Therefore, the 
precession is no longer around a common direction for 
all modes. If we started with a situation of maximum 
mixing, then the medium reduces the effective mixing 
angle for all modes. For a very dense medium, the ef- 
fective magnetic field will be almost perfectly along the 
2-direction, suppressing flavor oscillations entirely. 

This is very different if we consider a neutrino ensem- 
ble so dense that the neutrinos themselves produce a 
significant refractive index. In that case the medium's 
contribution to the refractive index is not along the fia- 
vor direction (i.e. along the z-axis), but rather along the 
direction of J. Put another way, neutrinos produce an 
"off-diagonal refractive index" because a given back- 
ground neutrino may be a coherent superposition of fla- 
vor states. The equation of motion for a single mode j 
now reads lllq 



B X Pj- + J X P 



2pj 



V 



J' 



(5) 



where the second term represents the self-interactions. If 
we sum this equation over all modes, then the second 
term becomes proportional to J x J = so that 



A 2 

J = B X > P,. 

^ 2p, ' 



(6) 



Therefore, the first derivative J of the ensemble's po- 
larization vector is not affected by the self-interactions. 
Still, the evolution of J is changed because the evolution 
of the individual modes is affected. However, the vacuum 
oscillations are not obviously suppressed even in a dense 
gas of neutrinos, in contrast to a standard background 
medium. 

Our case of active-active neutrino oscillations is very 
different from the active-sterile case Sterile neu- 

trinos do not produce a weak potential so that there is no 
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off-diagonal refractive index. The neutrino contribution 
to the "effective magnetic field" is along the ^-direction, 
i.e. not proportional to J. Therefore, the self-interaction 
term is somewhat more similar to the effect of an external 
background medium, although the oscillation equations, 
of course, remain non-linear. 



III. EXPLANATION OF SYNCHRONIZED 
OSCILLATIONS 

It is now easy to demonstrate that Eq. implies a 
synchronized precession of all modes around B if the neu- 
trino density is sufficiently large. To this end we first 
imagine the vacuum oscillation term to be absent, i.e. we 
consider the equation of motion 



P. 



J X P,. 



(7) 



Every individual mode precesses around the direction of 
J. Of course, if all neutrinos were initially prepared in a 
specific flavor state, then all Pj as well as J are aligned 
along the z-direction, and no precession takes place (per- 
fectly coherent state). Likewise, if the individual Pj 
point in random directions so that J = 0, again there 
are no precessions (perfectly incoherent flavor mixture). 
We consider the general case where the Pj initially point 
in many different directions, but do not add to zero. 

Next we switch on the vacuum term from Eq. (||), i.e. 
a weak external "magnetic field" . With "weak" we mean 
that for a typical mode the precession frequency around 
J is much larger than the one around B. This implies 
that the evolution of a given mode remains dominated by 
J. The fast precession around J implies that the projec- 
tion of Pj on J is conserved while the transverse compo- 
nent averages to zero on a fast time scale relative to the 
slow precession around B. If J moves slowly, then the 
individual modes will follow J. Put another way, the in- 
dividual modes are coupled to each other by their strong 
"internal magnetic fields," forming a compound system 
with one large magnetic moment. It is this compound 
object which precesses around B. Of course, if the ex- 
ternal field is much larger than the internal ones (dilute 
neutrino gas), then the modes will decouple and precess 
individually around the external field with their separate 
vacuum oscillation frequencies. 

Our picture is perfectly analogous to the Zeeman effect 
in atoms. An atomic state is characterized by its spin 
angular momentum S, its orbital angular momentum L, 
and the total angular momentum J = L -I- S. In a weak 
external magnetic field the spin-orbit coupling caused by 
the internal magnetic fields (Russell-Saunders coupling) 
remains intact, the external field is only a perturbation. 
In this case it is the total angular momentum J which 
precesses, i.e. which determines the atomic level splittings 
caused by the external S-field. On the other hand, if the 



external field is much stronger than the internal one, then 
L and S decouple and precess independently around the 
external field: the atomic levels are determined by the 
separate quantization of L and S along the external in- 
field (Paschen-Back effect). 

Granting that in the neutrino case the polarization vec- 
tors of the individual modes are locked in the sense that J 
indeed forms one large "angular momentum," the evolu- 
tion of the compound system is governed by the equation 



'-'synch 



B X J. 



(8) 



It remains to determine the value of CLigynch which plays 
the role of the gyromagnetic ratio for the compound 
system. Of the individual modes, the external field 
"sees" only the projection along J because the trans- 
verse components average to zero. Therefore, the con- 
tribution of mode j to the total magnetic moment M is 
{Am^/2pj)3 



Pj so that 



(9) 



where J is a unit vector in the direction of J. Therefore, 
the gyromagnetic ratio is 



IMI 



1 



•^synch 



|J| |J|^ 2p, 



^^^J.P, 



(10) 



In particular, if all modes started aligned (coherent flavor 
state) then |J| — N^, and J • Pj = 1 so that 



-"synch 



Ar 



2p 



1 



3 = i 



2p, 



(11) 



in agreement with 

If an external medium is present, the uncoupled modes 
precess around different B-field vectors rather than a 
common direction. When the modes are coupled by self- 
interactions, it is still the one J which precesses around 
one common B field which is a suitable average of the 
individual Bj which easily can be worked out. 

Returning to the simpler case of a common B for all 
modes, the calculation of Wsynch amounts to determin- 
ing the Lande-factor for J in the atomic analogy. This 
is the problem of calculating the magnetic moment of a 
system if the angular momentum is the sum of individual 
components which have magnetic moments with different 
gyromagnetic ratios. In this case the vector sums of the 
angular momenta and of the magnetic moments are not 
co-linear. In atomic physics, the spin angular momen- 
tum produces twice the magnetic moment of the orbital 
angular momentum, hence the complication. 

We stress that, contrary to the previous literature, our 
analysis shows that there is nothing special about the ini- 
tially aligned state, even though this state may be most 
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FIG. 2. Evolution of the Ve survival probability for three 
values of the neutrino momenta in the presence of a strong 
neutrino self-potential term (k = 10). 



motivated by the neutrino applications when all of them 
start in one flavor state. Any initial configuration of J 
precesses as one vector. Again, in atomic physics L and 
S can be combined in different ways to form one J. For 
example, a p-state with L = \ and S" = 1/2 can combine 
toaJ = 3/2oraJ=l/2 state; there is nothing special 
about the "aligned" state (J = 3/2). 

Moreover, even though the initial J precesses as a com- 
pound system, the individual modes do not oscillate in 
unison, except for the special case of perfect initial align- 
ment and infinitely strong self-coupling. In the general 
case the motion of every is a fast precession around 
J, superimposed on a slow precession of J around B. 
The compound motion of Pj is generally rather compli- 
cated and different for every mode. This is illustrated in 
Fig. H, where the synchronized oscillations of three differ- 
ent modes are shown for k = 10, starting with a perfectly 
ahgned state. 

Our treatment predicts that any set of initial Pj leads 
to synchronized oscillations in the sense of a precession of 
the initial J with a single frequency Wgynch- The length 
of the initial-state J is conserved while the oscillation 
frequency depends on details of the initial distribution 
of polarization vectors. The system does not prefer one 
particular synchronized state over another. For example, 
the completely incoherent state will stay that way, and 
will not spontaneously align the Pj to form a coherent 
state, in agreement with the stability analysis of ||]. Any 
initially prepared J will do nothing but precess about B. 
Put another way, the nonlinear aspect of the neutrino 
system manifests itself in the coupling of the individual 
Pj to each other to form a compound J which acts as 
one large angular momentum. Beyond this, the system is 
easily understood in terms of a linear equation of motion. 

Our analysis is also entirely independent of the num- 
ber of neutrinos or modes. For example, if there are only 
two modes nothing in our analysis changes so that the 
synchronization effect should not be viewed as a collec- 
tive phenomenon. In fact, the atomic example was one 



consisting of two coupled magnetic moments, the orbital 
and spin terms. We have checked with a simple numer- 
ical code that a system of two or three modes indeed 
behaves as expected according to our treatment. Actu- 
ally, a two-mode system can be fully solved analytically 
so that in principle the physical essence of the synchro- 
nization effect can be calculated without recourse to nu- 
merical methods. 

There are exceptions to our general statements which 
will be of interest in the next section. Evidently we can 
construct "pathological" cases where several angular mo- 
menta Pj add up to a vanishing or arbitrarily small J 
while the magnetic moments do not, leading to a gy- 
romagnetic ratio which can be constructed to become 
arbitrarily large. In this case Wgynch will no longer rep- 
resent some typical Ato^/2p, but can be constructed to 
become arbitrarily large. In this case the motion of J is 
not necessarily slow compared to the internal precessions 
of the individual Pj so that our treatment is no longer 
adequate. However, the perfectly incoherent state of an 
ensemble of many modes will not be pathological in this 
sense because the random distribution of P,- will ensure 



that both |J| - |M| - l/^ 
of modes. In the limit N^, 
simultaneously. 



N^, where Nv is the number 
CO both J and M vanish 



IV. NEUTRINOS PLUS ANTINEUTRINOS 

As a next step we may extend our analysis to the 
case where neutrinos and antineutrinos are simultane- 
ously present, an inevitable situation in a realistic system 
such as the early universe unless the neutrino chemical 
potentials are extremely large. To first order in Gf the 
equations of motion are 



P. 



-'bxP, + ^(J-J)xP, 



2W 



V 



™'bxP, + 4?^(J-J)xP,, 



2pfe 



V 



(12) 



where overbarred quantities refer to anti- neutrinos. In 
particular, the total polarization vectors are 



and 



fe=l 



(13) 



As explained in | p6[ , the definition of the polarization 
vector for anti-neutrinos is "reversed" in the sense that 
in vacuum it precesses in the opposite direction of neu- 
trinos This corresponds to fcrmions with a true mag- 
netic moment which have opposite gyromagnetic ratios 



*In the flavor oscillation case one can avoid this minus sign 
by defining both polarization vectors in the same way so that 
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for particles and anti-particles. For example, a neutron 
and an anti-neutron spin-precess in opposite directions 
in the same external magnetic field. 

It is obvious that a system consisting of neutrinos and 
anti-neutrinos behaves the same way as one consisting of 
neutrinos only, except that the role of the total angular 
momentum is now played by I = J — J. The anti-particles 
appear as normal modes of the system, except that they 
sport negative gyromagnetic ratios. It is I = J — J which 
precesses slowly around the external field, while all Pj 
and Pfe remain pinned to I. The corresponding evolution 
equation for the compound system is 



'^synch 



B X I. 



where the total gyromagnetic ratio is 



^synch 




I P. 



fc=l 



2pk 



(14) 



(15) 



We have checked with a numerical code several situations 
with a thermal population of neutrinos and antineutri- 
nos, and the results were always as expected. 

Therefore, a system consisting of neutrinos and anti- 
neutrinos typically behaves qualitatively similar to the 
neutrino-only case. However, the negative gyromagnetic 
ratios of anti-neutrinos relative to neutrinos allow for 
"pathological" situations where the system behaves qual- 
itatively differently than above. 

One such case is a thermal ensemble without chemical 
potential, i.e. a situation where N^, = A^p. If all neutrinos 
start in a given flavor state we have I = J — J = 0. This 
situation is "pathological" in the sense that a vanishing or 
very small I is associated with a large magnetic moment 
because particles and anti-particles enter with exactly 
opposite magnetic moments. To illustrate this case we 
consider only one mode of particles and anti-particles so 
that the equations of motion are 



the vacuum oscillation equations look the same. In this case 
the sign change has to be introduced by hand in the non-linear 
potential terms by flipping the sign of one component of the 
potential. This can be achieved by defining P* as the vector 
with the second component reversed, leading to equations of 
motion of the form 



{3 -3 )xP„ 



Pft = +- B X Pk - -^-TT-^ (J' - J) X Pfc. 



2pfe 



V 



This approach was used in much of the literature. While it is 
mathematically equivalent to our treatment, we think that it 
obscures the simplicity of the equations because their vector 
form is destroyed. 
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FIG. 3. Evolution of and jl| for the system of Eqs. (0), 
when sin 29 — 0.01 and in the presence of a strong neutrino 
self-potential with lu — —0.01, corresponding to an inverted 
mass hierarchy Am^ < 0. 



P = +CJ B X P + (P - P) X P, 
F = -cj B X P + (P - P) X P, 
with a suitable u. This implies 

P F==wBx (P + P). 



(16) 



(17) 



Evidently I = P — P is not conserved if the two polar- 
ization vectors start aligned so that at first 1 = 0. The 
effect of the external field is to drive P and P apart, 
creating a net I which is orthogonal to B. In the case of 
large mixing both P and P will oscillate much faster than 
they would in vacuum, yet convert to the other flavor. 

In the case of small mixing the result strongly depends 
on the sign of w in Eqs. (p^), i.e. on the sign of Am?. For 
LU > 0, the directions of the vectors P, P and B almost 
coincide. The role of the strong neutrino self-potential 
term is just to increase the oscillation frequency, while 
the amplitude of Pz is the same as in the vacuum case. 

For an inverted mass hierarchy (Am^ < 0) and small 
mixing angle, B is close to the z-axis, but it is almost op- 
posite to the initial directions of P and P. A small seed 
of |I| ^ is enough to drive P and P to the opposite di- 
rection from their initial orientation, i.e. one can achieve 
complete flavor conversion. We illustrate this situation 
in Fig. 1^, where the z-component of the polarization vec- 
tor and the modulus of I are plotted for w = —0.01 and 
sin 20 = 0.01. One can see that Pz evolves from 1 to 
— 1 and back. This resonance has a different origin from 
the Mikheev-Smirnov-Wolfenstein mechanism because it 
comes from the neutrino self-potential term. 

These phenomena were discovered and discussed in 
and subsequent papers of that series. We merely stress 
that our way of writing the equations allows for a 
straightforward visualization of the motion of the polar- 
ization vectors. 

Another physically motivated "pathological" case is 
when two flavors of neutrinos are thermally populated 
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with large but opposite chemical potentials. In the early 
universe this corresponds to a hypothetical initial condi- 
tion of a vanishing or small lepton-number density, yet 
an anomalously large density of radiation in the form of 
neutrinos. Since the chemical potentials are assumed to 
be opposite we have J = —J so that I = J — J = 2J 
is now "large." On the other hand Wgynch = so that 
flavor oscillations take place with a vanishing frequency, 
i.e. the initial condition is frozen without further evolu- 
tion. This appears to be the only case where the large 
neutrino-neutrino self-potential acts to prevent flavor os- 
cillations. 

These "pathological" cases have an atomic counterpart 
in the form of positronium where the two spins are associ- 
ated with equal but opposite magnetic moments. There- 
fore, the 1=1 state (ortho-positronium) has a vanishing 
magnetic moment and thus shows no Zeeman splittings 
in a weak external field. Likewise, the 7 = state (para- 
positronium) consists of only one level and thus cannot 
split. Therefore, even though we have a system of two 
spins associated with two magnetic moments, there is no 
weak-field Zeeman effect. In a strong field the electron 
and positron spins and magnetic moments are separately 
quantized along the external i3-field, giving rise to a non- 
trivial level structure. 



V. CONCLUSIONS 

In summary, we have provided a simple and physi- 
cal explanation of the synchronized oscillations observed 
in the numerical treatment of dense neutrino ensembles. 
The effect is perfectly analogous to the coupling of sev- 
eral angular momenta, for example spin and orbital angu- 
lar momentum in an atom, to form one large compound 
angular momentum with one large associated magnetic 
dipole moment which precesses as one object in a weak 
external field. Anti-neutrinos are naturally included in 
our picture if one observes that they carry "negative gyro- 
magnetic ratios" in flavor space. The nonlinear nature of 
neutrinos oscillating in a background of neutrinos is thus 
reduced to a very simple and well-known coupling effect 
of magnetic moments to each other. Our picture provides 
a transparent framework that accounts perfectly for all 
of the previously discussed synchronization phenomena 
in the literature, and that allows one to derive both gen- 
eral properties and special cases in a unified framework. 

The practical impact of synchronized neutrino oscilla- 
tions in the early universe will be studied in a forthcom- 
ing paper 1 18|. 
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